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The ofF-resonant hyperpolarizability is calculated using the dipole-free sum-over-stats expression 
from a randomly chosen set of energies and transition dipole moments that are forced to be consistent 
with the sum rules. The process is repeated so that the distribution of hyperpolarizabilities can 
be determined. We find this distribution to be a cycloid-like function. In contrast to variational 
techniques that when applied to the potential energy function yield an intrinsic hyperpolarizability 
less than 0.71, our Monte Carlo method yields values that approach unity. While many transition 
dipole moments are large when the calculated hyperpolarizability is near the fundamental limit, 
only two excited states dominate the hyperpolarizability - consistent with the three-level ansatz. 



PACS numbers: 42.65, 33.15.K, 33.55, 42.65.A 



I. INTRODUCTION 

Large nonlinear-optical susceptibility materials 
are central for many optical applications such 
as telecommunications, [H three-dimensional nano- 
photolithography,0,Q and new materials development [1| 
for novel cancer therapies. 01 In the present work, we 
focus on the molecular hyperpolarizability, which is the 
microscopic basis of electro-optic switches and frequency 
doublers. 

The fundamental limits of the hyperpolarizability had 
been calculated using what appeared to be three-level 
truncation of the sum rules. [1, Q Champagne and Kirt- 
man were critical of the accuracy of truncated sum 
rules. In response, the three- level ansatz was intro- 
duced, which states that when a quantum system has 
a hyperpolarizability near the fundamental limit, the 
system can be described by a three-level model. Q Us- 
ing a new dipole-free form of the sum-over-states ex- 
pression for the hyperpolarizability, [l3| Kuzyk showed 
rigourously that the three-level ansatz does not demand 
the sum rules to be truncated. TT| Furthermore, the 
three-level ansatz with the new dipole free expression was 
used to calculate the more general fundamental limits 
of the dispersion of the real and imaginary parts of the 
hyperpolarizability. [UJ This result can be used to study 
any phenomena that is based on a second-order response. 
More importantly, since the fundamental limit can be 
determined at any wavelength, it can be used as an ab- 
solute standard for comparing hyperpolarizabilities. In 
particular, the ratio of the measured hyperpolarizability 
to the fundamental limit, called the intrinsic hyperpolar- 
izability, can be used to compare the efficacy of any two 
molecules, independent of their size. 

A factor of thirty gap between the best molecules 
and the fundamental limit has persisted even after three 
decades of molecular engineeriiigjTj and there is no 
known explanation for this gap. [12] It is not of a fun- 
damental nature: a simple clipped harmonic potential 
energy function, which can be solved analytically, falls 
well into the gap. As such, a numerical optimization 
technique, which varies the potential energy function in 



a way that maximizes the intrinsic hyperpolarizability, 
was introduced to yield a better understanding of the 
nature of what makes a large hyperpolarizability. 13] In 
that work, Zhou and coworkers found that an oscillat- 
ing wavefunction serves to localize the wavefunctions in 
a way that forces the system to have only two dominant 
excited states that interact with the ground state, jisj 
This result showed how the three-level ansatz is a nat- 
ural way for optimizing the intrinsic hyperpolarizability. 
Based on this observation, Zhou and coworkers proposed 
that molecules with modulation of conjugation, i.e. with 
wiggly wavefunctions, may be a new paradigm for making 
molecules with exceptionally large hyperpolarizabilities. 

Perez Moreno and coworkers reported on a series of 
molecules with varying degrees of modulation in the con- 
jugated bridge between the donor and acceptor ends 
of these quasi-one-dimensional molecules. |U] Molecules 
with uneven aromaticity in the bridge, as suggested by 
Zhou's work, were measured to have larger intrinsic hy- 
perpolarizability than those with a smoother bridge. The 
best molecule was reported to have a world's record in- 
trinsic hyperpolarizability - 50% better than the previous 
best. 

Zhou and coworkers continued to study various ways 
of optimizing the potential energy function and found 
that modulated conjugation was but one way of getting 
near the limit. jTsj Depending on the starting potential 
given to the optimization program, the optimized poten- 
tial can be wiggly or smooth; the wavefunctions for the 
optimized potential can be strongly overlapping or well 
separated; and only a few states or many can contribute. 
In all cases, however, the numerical value of the intrinsic 
hyperpolarizability is just below 0.71 - about 30% below 
the fundamental limit. Thus, there appears to be a large 
number of distinct potential energy functions that yield 
the same maximum hyperpolarizability. 

These results clearly show that while the three-level 
ansatz appears to accurately the upper bound of the hy- 
perpolarizability, systems that are not well approximated 
by three-levels can also have an equally large nonlinear 
response. But, it is puzzling that all calculated locally 
maximal values fall exactly 30% short of the limit. In 



2 



the numerical optimization studies, the potential energy 
function is varied and the transition moments and ener- 
gies calculated from the resulting wavefunctions. Could 
the sum rules allow for transition moments and energies 
that do not result from such simple potential energy func- 
tions? 

To answer these questions, we take the approach of 
varying the transition moments and energies directly 
while enforcing compliance with the sum rules, circum- 
venting the need for calculating wavefunctions from a 
potential. This process is numerically more efficient and 
faster than solving the eigenvalue problem for a Hamil- 
tonian, and allows us to investigate a much larger pa- 
rameter space. The downside is that we cannot easily 
associate a potential energy function with those matrix 
elements and energies that yield the largest hyperpolar- 
izability. 



II. APPROACH 



We express the sum rules in dimensionless form, 
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(1) 



where e„ = Eno/Eio and ^np — Xnp/ x^j/^-^ . En is the 
energy of state n and E^m = — E^- —eXnm is the 
electric transition dipole moment between state n and 



and 
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h Ne/2mEiQ is the fundamental limit 



of the transition moment between the ground and first 
excited state, where m is the mass of the electron, e is 
the magnitude of the electron charge and iVg the number 
of electrons. Equation [T] represents an infinite number of 
equations, one for each distinct value of m and p. We 
can therefore refer to each distinct equation by the pair 
of indices {m,p). 

In the limit of an [N + l)-level model (i.e. only N 
excited states contribute to the SOS expression for the 
hyperpolarizability), we can treat (ei, . . . , e^v) and 
for n, m < N as parameters. Note that we label the 
ground state with the subscript zero. We pick the values 
of these parameters randomly, but constrained by the 
sum rules, as follows. 

We start by defining the energy levels. For an (A^ + 1)- 
level system, we pick N random numbers in the range 
< r < 1 and order these in a vector in ascending order 
so that ri < r2 • ■ • < rjv. It is important that ri 7^ since 
all of the dimensionless energies are given by e„ = rn/ri. 
Note that cq = and ei = 1. 

Once the energies are determined, the sum rules given 
by Equation [T] are used to fix the transition moments, 
as follows. Beginning with the sum rule Equation (0,0), 
we start by picking a random number, r, in the range 
— 1 < r < +1. Since each term in the sum in Equation 
(0,0) is non-negative, l^oil^ei < 1- But, since ei — 1, 
this implies that |^oi| < 1- We assume that all of the 
matrix elements of the dipole operator are real, so we set 



foi = ^10 = f- Equation [T] with N excited states and 
m = p = Q can be written as: 
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^ (e„ - eo) C™o = 1 - (ei - eo) i 
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n=2 



where the numerical value of (ei — eo) Cio is computed 
from the already-determined values of ei and fio. Note 
that Co = 0, so we drop it. Since Equation [2] implies that 



^lo < (1 - ei^i'o) /e2, 



(3) 



we pick a random number r' in the range — 1 < r' < +1 
and get ^02 using the expression, 



^02 = 60 = J''\/(l-eiCio)/e2- 
To get ^30, we evaluate 



Co3 = 60 = ''"\/(l-eiC?o-e2^|o)/e3, 



(4) 



(5) 



were r" is another random number. Thus, for the set 
of random numbers r^r' ^r" ^ . . we continue the process 
until we have the matrix element ^10, 'C201 Csoj • ■ • Cwo- 

Next, we use the sum rule equation (1,1). Using the 
same procedure as above, for a random number — 1 < 
s < +1, 



62 = 61 = sV^(l + eiC?o)/(e2-ei). (6) 
For -1 < s' < +1, 

Ci3 = 61 = s'^ (1 + ei^io - (62 - ei)Cii) /(ea - ei), 

(7) 

and so on, until we get £,ni- We continue the process for 
sum rules (n, n) , for all n where n < N. This procedure 
then yields all of the off-diagonal transition moments. 
Note that we ignore the sum rule {N, N) because it gives 
the nonsensical result that the oscillator strength is neg- 
ative. 

With this protocol, transition moments of the form ^nm 
for m ^ n and the energy levels are determined by the 
role of dice; but, these parameters are forced to obey the 
sum rules. We then use these transition moments and 
energy levels to calculate the dipole-free intrinsic hyper- 
polarizability. Pint, which is given by:pi)| 



Pint — 



P 



P 



MAX 



3\ 
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2e, 



(8) 

where Pmax is the fundamental limit of the hyperpolar- 
izability and P is the hyperpolarizability. 

Note that because the dipole-free expression was de- 
rived from the sum rules, and explicitly eliminates dipole 
moments of the form ^„„, we do not need to determine 
electric dipole moments of the system's eigenstates to 
calculate the intrinsic hyperpolarizability. 
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III. RESULTS AND DISCUSSION 

Past studies of optimizing the hyperpolarizabilities 
used the potential energy function as a starting point. (Tsl. 
[l5j The wavefunctions for a particular potential energy 
function are determined from the Schrodinger Equation. 
These are then used to determine the dipole matrix and 
energies of the system, which are in turn the input pa- 
rameters for calculating the hyperpolarizability using the 
sum-over-states expression. 

The potential energy approach yields several avenues 
for investigating properties of a system that yield the 
largest hyperpolarizability. For example, one can arbi- 
trarily vary the potential energy function to determine 
directly how the shape of a molecule, characterized by the 
potential, affects the optical response. Alternatively, one 
can define a coulomb potential energy function for var- 
ious arrangements of nuclei to study structure-property 
relationships, [l^ The latter approach, while more real- 
istic because it more closely describes real molecules, is 
more limited since it is restricted to a superposition of 
a finite number of point sources. The former case, on 
the other hand, might lead to prescriptions for making 
artificial systems such as multiple quantum well layered 
structures or nano-engineered structures. 

A more basic question pertains to the limitations im- 
posed by restricting the Hamiltonian to be of the form 
of a kinetic energy term and a potential energy func- 
tion that depends on position. It is possible that certain 
combinations of energy levels and dipole matrices may 
be derivable only from a more general Hamiltonian that, 
for example, contains terms such as p - f{r) + f{r) -p (i.e. 
products of functions of the position, r, and the momen- 
tum p) but that may never-the-less still obey the sum 
rules. Perhaps these more general Hamiltonians fall into 
classes that each yield a different fundamental limit, but 
with Pint < 1. Our Monte Carlo method does not dis- 
criminate between these different Hamiltonians as long as 
the sum rules are obeyed. Alternatively, matrix elements 
and energies that are consistent with the sum rules may 
be inconsistent with some additional basic physical prin- 
ciple. These are interesting questions that will not be 
explored here. Rather, we focus on investigating what 
matrix elements and energy levels lead to an improved 
nonlinear optical response; and, how these results com- 
pare with those obtained from a simple potential energy 
function, V{r). 

Figure [IJi shows the distribution of the calculated in- 
trinsic hyperpolarizabilities for 10,000 runs using a 5-level 
model. First we note that all of the calculated values of 
the intrinsic hyperpolarizability have a magnitude that 
is less than unity. This is consistent with the existence 
of a fundamental limit. Furthermore, the largest values 
of Pint approach unity. In contrast, when the potential 
energy function is varied, [ll] independent of the starting 
function, the maximized values converge to about 0.7. 
So, our results suggest that the hyperpolarizabilities at- 
tainable from the simple potential energy functions used 
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FIG. 1: Distribution of calculated hyperpolarizabilities for 
various energy weightings. Error bars are calculated from the 
square root of the frequency of observing a particular value 
of Pint- 



in past studies may not be as large as is possible based 
on only the sum rules. How the sum rules, which are 
derived from commutation relationships that assume the 
potential energy function to be dependent on only the 
positions of the electrons, yield a different upper bound 
than direct use of the potential energy function is an open 
question that will be the focus of future studies. Perhaps 
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FIG. 2: The matrix elements that yield the largest value of 
the intrinsic hyperpolarizability {Pi„t ~ 0.9569) after 100,000 
tries for a 10-state model. 




FIG. 3: The fractional contribution, /9^™, for a 10-state 
model with Pint ~ 0.9569 (using matrix elements from Figure 

ID). 



such investigations will shed light on why all molecules 
fall a factor of 10'^/^ below the fundamental limit .[l^ 

If the energy levels are chosen randomly, they will on 
average be equally spaced. To probe a larger parameter 
space, we define a power law weighting function when 
randomly choosing the energy level spacing. Thus, on 
average, the energy level spacing will be of the form 
Figure [T^ shows the result for weighting factors with ex- 
ponents p = 2 and p = 0.5. The solid curve shows the 
function 



F = A 



Pint) 



(9) 



which is fit to the p — 2 simulation with n and A as 
adjustable parameters. The best fit yields n = 2. 1.38% 
of the values fall outside the range |/3mt| > 0.71. 

Figures [T)d and [T]; show the 5-level model results for 
p = 0.5, 0.75, 1.0, 1.5, and 2.0. The curve given by Equa- 
tion [5] with n = 2 is shown for reference. Clearly, the 
simulations with smaller weighting factors, p, yield a nar- 
rower distribution, showing that the number of simula- 
tions giving large intrinsic hyperpolarizability is smaller 
then the case when p is larger. This is consistent with 
our recent 3-level model analysis that suggests the sec- 
ond excited state energy should be much higher than the 
first in order to optimize the hyperpolarizability, [l^, [13 
that is, small E10/E20 leads to large hyperpolarizabil- 
ity. Since, on average, E10/E2Q = (1/2)^ - the intrinsic 
hyperpolarizability is larger for larger p. 

To study the character of the states of a system 
that leads to the largest nonlinear response, the tran- 
sition dipole moments and energies that yield the largest 
value of the intrinsic hyperpolarizability from the set of 
100,000 random values generated by our protocol is stud- 
ied in more detail. Figure [2] shows the matrix elements 
corresponding to the largest hyperpolarizability(/3i„i = 
0.9569) in the set for a ten-state model. For this set. 



the energies, in dimensionless units of the form e„ = 
{E,, - Eo)/Ei are : 0, 1, 10.08418, 11.16255, 17.29753, 
20.40493, 22.37591, 27.98258, 34.93178, and 45.98839. 

There are several excited states with substantial tran- 
sition moment, though most of the oscillator strength 
resides in the lower-lying states. However, in assessing 
the contribution of each state to the hyperpolarizability, 
it is more appropriate to use Equation [8] to define the 
fractional contribution from pairs of individual states n 
and m, 



Pint ~ 1 



3/4 



1 



2e, 



(10) 



Figure [3] shows a plot of the fractional contribution to 
the intrinsic hyperpolarizability due to states n and m. 
Note that which is calculated from Equation [TUl is 

not symmetric to interchange of n and m. However, since 
(3 is calculated from a sum over terms of the form (3^^ + 
/3™j", the resulting hyperpolarizability is symmetric upon 
interchange of the indices. 

Clearly, the two states that contribute most strongly 
to the intrinsic hyperpolarizability are states 1 and 2. 
State 3 contributes, but is small in comparison. All other 
states are negligible. This clearly shows that when the 
hyperpolarizability is maximized, the system reduces to a 
three-state model - consistent with the three-level ansatz. 

It is an interesting exercise to assess this large- 
hyperpolarizability "molecule" using an analysis pro- 
posed to treat complex molecules in terms of parameters 
determined from properties of the first two dominant ex- 
cited states. [12] In this analysis, the sum rules reduce the 
SOS expression in the three-level ansatz top^ 



Pint = f{E)G{X), 
f{E) = (1 - Ef^ 



where E = Ew/E^q, X = xiq/.t^^^^ 



E 



(11) 
(12) 
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FIG. 4: (a) The matrix elements that yield the largest value of 
the intrinsic hyperpolarizability {Pi„t ~ 0.9441) after 100,000 
tries for a 40-state model; and (b) the fractional contribution, 

n,m 
int ' 

and 

GiX) = m^^il-X'^). (13) 

Equation [11] should be a good approximation to a 
molecule whose nonlinear response is close to the limit. 

With E = 1/10.08, we get /(0.099) = 0.991. The 
transition moment to the first excited state is given by 
X = 0.712, which yields G(-0.712) = 0.989. Thus, 
the intrinsic hyperpolarizability \Pint\ = 0.98, which is 
in excellent agreement with /3int = —0.96 when all 10 
states are used to calculate the intrinsic hyperpolarizabil- 
ity. Note that while many transition moments appear to 
be substantial (as we can see from Figure ^ , the three- 
level ansatz yields a good approximation. 

The same pattern holds when more states are included. 
Figure shows the matrix elements of a 40-state model 
that yields Pint = 0.9441. This is the largest value after 
100,000 trials. As in the 10-level case, we see lots of 
large matrix elements; but, the fractional contribution 



fHiNT Dominant States Number of Levels 

0.033 0, 1, 2, 3, 5, 7, and 8 7 

0.31 0, 1, 3, 4, and 7 5 

0.95 0, 1, and 3 3 

TABLE I: States whose contributions are at least 10% of the 
largest contribution for the 40-state model shown in Figure |4l 

to the intrinsic hyperpolarizability, as shown in Figure 
[IJa is dominated by one term corresponding to the first 
and third excited states. Again, the three-level ansatz 
appears to hold. 

The relevant state are characterized by A = 0.799 and 
E = 0.0433, which yields (3int = /(0.0433)G(0.799) = 
0.998 X 0.991 = 0.99. This is in reasonable agreement 
with f3iNT = 0.944, which is obtained when all states are 
included in the sum. In contrast to the 10-level model, 
the 40-level model has (40/10)^ = 16 times more con- 
tributions to the hyperpolarizability. So, when all terms 
are included there is a larger chance for many smallers 
terms to add to a significant amount. 

In all of our calculations, when the number of states are 
increased, for example, to 80 states and beyond, the same 
pattern appears to hold: while many matrix elements 
are large, only two excited states dominate the intrinsic 
hyperpolarizability. 

Figure [5] shows the contributions to the hyperpolar- 
izability for two systems, with hyperpolarizabilities of 
0.033 and 0.31. More states contribute as the intrinsic 
hyperpolarizability gets smaller. Note that for these two 
cases, similar to the above results, many transition mo- 
ments are large, but only a few combinations contribute 
to (3. 

To quantify this observation, we define a term in the 
SOS expression for P^j^ to contribute substantially if its 
magnitude is at least 10% of the term with the largest 
contribution. Table [J summarizes the dominant states 
that correspond to the terms that contribute substan- 
tially to the hyperpolarizability. Indeed, the number of 
states that contribute is typically larger when the intrin- 
sic hyperpolarizability is lower. We stress that this ob- 
servation is statistically based. There are certainly cases 
where the hyperpolarizability can be large even when 
many states contribute, as we have found in optimiza- 
tion studies of the potential energy function. [15] How- 
ever, we conclude that when we randomly pick transition 
moments and energies, the intrinsic hyperpolarizability is 
usually small when many states contribute. 

It is trivial to show from the sum rules that contribu- 
tions to /3 from states of higher energies can dominate 
over the lower energy states. However, since the three- 
level ansatz implies that the higher-energy states are not 
important when the hyperpolarizability is near the fun- 
damental limit, the higher- lying states must always can- 
cel each other if any of the terms are large. 

In our protocol, we force the system to be an N-level 
model. We can study the effects of the truncation of 
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FIG. 5: The fractional contribution, for a 40-state 

model when the intrinsic hyperpolarizability is 0.033 (top) 
and 0.31 (bottom). 

the sum rules by considering how the hyperpolarizability 
changes when we add an additional state. Recall that 
the diagonal sum rules when truncated to N-levels are of 
the form, 

N-l 

Tt = 

where we use the shorthand notation e„p = e„ — ep. 
Equation [HI holds for all p < — 1, where we find that 
p = N — 1 yields a nonsensical result, so that equation is 
ignored. 

If we add one more state, the diagonal sum rules for 
p < N — 1 are of the form, 

N 
n=0 

With the help of Equation [TH Equation [15] reduces to 
^N,p^p,NeN,p = 0. (16) 



Assuming that the added state is not degenerate with 
any of the other N states, this yields 

^N,p = £,p,N = 0. (17) 

Next we consider the case with p = N ^ 1. The sum 
rule, 

N 

en,N-l^N-l,p^p,N-l = 1 (18) 

with the fact that en,N-i = —&N-i.n < reduces to, 

N-l 

eN,N-l^N,N-l^N-l,N = 1 + ^ eN~l,n^N-l,n^n,N-l 

ri=0 

(19) 

So, it must be that ^Ar,Ar_i and eM,N-i are non-zero and 
are related to the values of the transition moments and 
energies that were determined by the Monte Carlo as- 
signments. 

The terms in the dipole-free expression for the hyper- 
polarizability, given by Equation IIOI have numerators of 
the form ^oi£.ij^jo, where i ^ j. Thus, the only additional 
terms that contribute to the hyperpolarizability that in- 
clude the state N are of the form £,q.n~i^n-i,n£,n,o and 
^o.p£,p.N£,N.o, wherep < A^ — 1. Both of these terms vanish 
because £,no = according to Equation [T71 

If we continue to add more states, we will find that 
transition moments between these new states and states 
TV, -I- 1, etc. will be nonzero; but, the transition 
moment between the ground state and the newly-added 
states will vanish because the oscillator strength is de- 
pleted by the lowest N states. The net result is that the 
additional states do not contribute to the hyperpolariz- 
ability. Since we can continue to add states in this man- 
ner, ad infinitum, to construct a full set of non-truncated 
sum rules, using truncated sum rules in conjunction with 
our Monte Carlo method does not appear to lead to any 
pathologies when calculating /3. 

One might ask if it is physical to have a quantum sys- 
tem in which transition moments from the ground state 
to all states outside of an N-level subspace vanish. The 
harmonic oscillator is an example of a system where tran- 
sitions from the ground state to all excited states beyond 
the first vanish. However, it may be that our Monte 
Carlo simulations may yield transition moments and en- 
ergies that, while consistent with the sum rules, may not 
be derivable from a potential energy function. 

Champagne and Kirtman had criticized sum-rule trun- 
cation, arguing that the resulting limits may be in 
error. Q The later development of the dipole-free SOS ex- 
pression for calculating the hypoerpolarizabilityfio'l lead 
to a more rigorous argument of why the original calcu- 
lations, which appeared to truncate the sum rules, did 
not violate the higher- level sum rules. pH A more impor- 
tant concern was the use of the three-level ansatz, which 
states that the fundamental limits can be calculated by 
assuming that only two excited states contribute. Given 
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that the sum rules allow transition moments between 
excited states to be large, and potentially many states 
can contribute, one might expect that the calculations 
of the fundamental limit would underestimate the upper 
bounds. This is clearly not the case. While a system 
with only two dominant excited states yields the funda- 
mental limit when the ratio of the oscillator strengths 
to each of the two excited states is chosen appropriately, 
systems in which many excited states contribute can also 
reach the fundamental limit. The implication is that all 
of the large contributions from transitions between ex- 
cited states cancel. 



IV. CONCLUSIONS 

We have shown that when assigning transition moment 
matrix elements and energies randomly, but under the 
constraint of the sum rules, values of the intrinsic hy- 
perpolarizability approach unity. This is in contrast to 
the intrinsic hyperpolarizability, which when calculated 
from an optimized potential energy function, falls bellow 
0.71. The distribution of the values of (iiNT is found to 
be well approximated by a cycloid-like function. When 
the energy-level spacing is weighted in a way that forces 
the lower energy states to have more closely-spaced en- 
ergies, the distribution falls below the cycloid function, 
suggesting that more sparse energy-level spacing of the 
lower levels yields a larger intrinsic hyperpolarizability - 
an observation consistent with past analytical studies. 

As the intrinsic hyperpolarizability increases, it is 
found that fewer states contribute even though the transi- 
tion moments between many states are nonzero and can 
be large. Furthermore, when /3/7vt is near unity, only 
three states contribute significantly, and the energy level 
spacing between the first and second excited states is 



large compared with the energy difference between the 
first excited state and the ground state. This result 
is consistent with the three- level ansatz. Furthermore, 
these matrix elements and energies yield an energy func- 
tion f{E) and transition moment function G{X) that ac- 
curately predicts the hyperpolarizability, which confirms 
the validity of the type of analysis proposed earlier, [l^ 
and applied to an amy lose helix. [17] 

We have also shown that our approach of using a fi- 
nite number of states, N , does not lead to pathologies. 
We did so by constructing a state space of dimension 
+ 1, whose transition moments and energies satisfy 
the higher-order sum rules; but, lead to the same hyper- 
polarizability as for the A^-dimensional state space. Since 
we can construct an infinite-dimensional space by succes- 
sively adding one state at a time, we conclude that the 
full set of sum rules are obeyed and the hyperpolarizabil- 
ity remains unchanged from the A^-dimensional one. 

We question whether or not the energy levels and tran- 
sition moments that are used to calculate the hyperpolar- 
izability when Pint > 0.71 are derivable from a simple 
Hamiltonian with a position-dependent potential energy 
function. It may very well be the case that, while the 
sum rules are derived from commutations relationships 
that assume a potential energy function that depends on 
position, using only the sum rules as a constraint to de- 
termining transition moments and energies may lead to 
results that belong to a class of more complex Hamil- 
tonians that also satisfy the commutation relationship 
[[if, a;], a;] = jm. If this is so, additional constraints 
may be needed to ensure that such systems are physi- 
cally reasonable. 
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